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Abstract 

Semileptonic, B — > ir W, decays are studied on quenched anisotropic lattices 
using tree-level tadpole improved Symanzik glue, NRQCD heavy quark and 
D234 light quark actions. Constrained fitting methods are applied to extract 
groundstate contributions to two-point and three-point correlators. We agree 
with previous lattice determinations of the form factors. The major source of 
systematic error here, as in previous work, comes from the chiral extrapolation 
to the physical pion mass. Future calculations must work at lighter quark 
masses to resolve this. 
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I. INTRODUCTION 



Determinations of the CKM matrix element \V u b\ from exclusive semileptonic B decays 
rely on lattice input for the (n\ J^B) or (p\ J^B) matrix elements. Reducing errors in these 
lattice calculations and extending their kinematic range will be an important contribution to 
consistency tests of the CKM matrix, and thereby of the Standard Model. Several quenched 
studies of B — > 7r IV decays have been carried out in recent years employing a wide range of 
different actions and methods fTl-ffl. Despite the differences in sources of systematic errors 
there is general agreement in final results, especially for the form factor f+{q 2 ) which is 
directly relevant for the differential decay rate measured in experiments. 

The difficulties with lattice simulations for these semileptonic form factors are well known. 
Lattice calculations are most reliable at large q 2 = {p^ —Pb) 2 , i-e. close to the zero recoil 
point, whereas experiments are limited to smaller q 2 . In the B meson rest frame, small q 2 
implies large pion momenta which introduces both ap n discretization errors on the lattice 
and large statistical errors. Furthermore it is becoming increasingly more evident that a 
major source of systematic error comes from the chiral extrapolation to a physical final 
state pion. 

This article describes further investigations of B — > 71 W decays in the quenched approxima- 
tion, where we have experimented with ways to improve errors due to finite pion momenta. 
The first difference from previous work is the use of more highly improved actions. We em- 
ploy Symanzik glue rather than the Wilson glue action, the D234 light quark action rather 
than clover and an NRQCD action corrected through 0(a 2 ). We find good dispersion re- 
lations for pion energies up to at least 1.5GeV. The ability or inability to cover a q 2 range 
overlapping with experiment is determined more by the other challenges, namely statistical 
noise and chiral extrapolations, than by discretization errors. This will be even more true 
if one goes to lattices finer than the relatively coarse, l/a s = 1.2GeV, lattices used in the 
present study. 

The other two new ingredients here are the use of an anisotropic lattice and constrained 
fitting methods. As demonstrated in many contexts, anisotropic lattices lead to correlators 
that have a larger number of data points inside the time region where signal-to-noise is still 
good. This allows for more accurate extraction of energies and amplitudes 0-0]. Reference 
II showed the advantages of anisotropic lattices for correlators involving finite momentum 
hadrons. In the present work we add one more major improvement to this scheme, namely 
constrained fits IJ. These new analysis tools, based on Bayesian statistics, allow us to 
increase the number of exponentials in fits to single correlators, without losing stability or 
having errors in low lying energies become large. One typically fits to all (or almost all) data 
points, with t min = or = 1, irrespective of where or whether a plateau sets in. Hence one 
can take full advantage of all the data points with small statistical errors that the anisotropic 
lattice gives us. Previously if one just went to anisotropic lattices and used conventional 
fitting methods one had to rely on having excellent smearings or a large matrix of smearings 
to ensure overlap between a plateau region and the region with good signal-to-noise. 

We find that constrained fits allow us to extract groundstate contributions to two-point 
and three-point correlators in a controlled way despite the fact that our smearings are far 
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from optimal. Nevertheless, our final results for B — > tx Iv form factors still exhibit large 
systematic errors, comparable to those quoted by other groups. The main reason for this is 
that our systematic error is dominated by the chiral extrapolation uncertainties. The D234 
light quark action suffers, in common with the clover and Wilson actions, from exceptional 
configurations. This limits our ability to go to small pion masses and necessitates a large 
chiral extrapolation (we work in the range 0.7m s t r ange < Tn q < ^-^f^strange)- Even the correct 
ansatz for the extrapolation is unclear at the present time. In the future, rather than try 
to go to finer lattices with the same action, we believe that, in order to improve on chiral 
extrapolation errors, one needs to go to light quark actions with good chiral properties. 



Work with improved staggered light quarks, for instance, has already started |[10|JTT|| . 
Difficulties with chiral extrapolations also limit our ability to go to pion energies larger than 
lGeV (or q 2 < 16GeV 2 ). If one were to work harder and double the statistics one might be 
able to consider pions with momenta up to 1.5GeV and go down to q 2 ~ 12Gev 2 as long 
as the light quark mass stays around the strange mass. Nevertheless, carrying out chiral 
extrapolations at large pion momenta would introduce prohibitively large errors. Recently 
the idea has been put forward of working in a reference frame where the entire q 2 range 



can be covered with pions having momenta less than lGeV (see ref. |12[ for discussion of 



the "Moving NRQCD" formalism). We believe the future of accurate semileptonic B decay 
studies lies in a combination of using light quark actions with good chiral properties and a 
formalism that can handle B mesons decaying at large velocities. The lessons derived from 
the present project should be very useful in such future work, especially in the analysis of 
three-point correlators. 

In the next section we introduce the actions used in this project. Section III provides 
simulation details and discusses results from two-point correlators. Section IV describes 
our direct fits to three-point correlators. This step replaces the conventional approach of 
considering ratios of three- and two-point correlators and looking for a plateau which is 
then fit to a constant. In section V we discuss chiral extrapolations and present results for 
form factors. Comparisons are made with previous lattice work. We then conclude with 
a summary section. We delegate details of constrained fitting methods to an appendix. A 
further appendix lists heavy-light current one-loop matching coefficients for the action used 
in this article. 



II. GAUGE AND QUARK ACTIONS 

The gauge and quark actions of this article are the anisotropic actions discussed in ref. [||. 
For the glue we use the tree-level tadpole-improved Symanzik action with rectangles only in 
spatial directions. 



1 J 5 P ss i 1 R ss i 1 R s / S 

,., s ,xTi3^r~i2^r~i2^r 




4 P st 1 Rst 



3 u 2 u 2 12 u 4 s u 2 t 
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P^ v and denote plaquettes and rectangles in the \w plane. The variables s and s' run 
only over spatial directions and u t and u s are the tadpole-improvement "wo" factors for 
temporal and spatial link variables respectively. We use the Landau link definition of Uq in 
this article. Xo is the bare anisotropy which differs from the true or renormalized anisotropy, 



X = a s /a t , 



(2) 



in the presence of quantum corrections. We use torelon dispersion relations |I3| to fix x 
nonperturbatively. Starting with xo — 3. at (3 = 2.4 we find the renormalized anisotropy to 
be x = 2-71(3). ' 



The light quark action is the D234 action of ref. |L4| modified for anisotropic lattices. 
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The quark fields \l/ c and the dimensionless lattice fields \1/ are related through 



(4) 



(5) 



Definitions of lattice derivatives and the improved are given, for instance, in the appendix 
to ref. |15| . We work with r, Cf, C3 and C4 all set equal to one. The "speed of light" 
coefficient, Co, is tuned perturbatively or nonperturbatively using pion dispersion relations. 
Its actual value will be discussed in the next section. 

For the heavy quark we use the standard NRQCD evolution equations which follow from 
the action IflMl, 



S 



NRQCD 



H is the nonrelativistic kinetic energy operator, 

A< 2 > 



a t 5H\ 



t-i 



$t-i 



a t H 



2 X (a s M y 

and S H includes 1/M relativistic and 0(a 2 ) finite lattice spacing corrections, 

1 - 



a t 5H 
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(6) 
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All derivatives are tadpole improved and, 

A (2) = Evf, A( 4 )=^vf (9) 

3=1 3=1 

The leading discretization errors in the total action are 0(a s a s ) errors coming from the 
light quark action. The leading relativistic corrections are 0(a s KQ ^ D ) coming from one- loop 
corrections to tree-level coefficients in the NRQCD action. We will work with heavy-light 
currents corrected to the same level as the action, i.e. we will include 0(a s ) and 0(Aqcd/M) 
terms but not 0(a s a s ) or 0(a s Kq ^ d ) corrections. We estimate systematic errors from the 
latter terms to be at the 8% and 3% levels respectively. 



III. SIMULATION DETAILS 

Table I summarizes lattice and action parameters. We work on 12 3 x 48 quenched anisotropic 
lattices with \ = a s/ a t — 2.71(3), as determined from torelon dispersion relations in the 
pure glue theory. We use a total of 199 configurations and run both time-forward and time- 
reversed NRQCD evolutions in order to increase statistics. The former uses timeslices - 
23, and the latter timeslices 0,47 ,46 ....25. For the B meson correlators and the three-point 
correlators we find little evidence for correlations between the two runs and close to a \[2 
improvement in statistics. Nevertheless, in our data analysis we always bin the data from 
the two time evolutions before carrying out fits. 

Based on string tension calculations we estimate the spatial lattice spacing to correspond to 
l/a s = 1.20(5)GeV. The p mass gives a similar value of (l/a s ) p = 1.18(6)GeV. 

We have carried out simulations at five values of the bare light quark mass corresponding 
to P/V = 0.624(13), 0.675(8), 0.714(14), 0.736(6) and 0.760(6). In terms of a t m = a t (m - 
m crit.), we have the range a t m = 0.023, 0.028, 0.033, 0.038 and 0.043. The middle value 
of a t m = 0.033 is very close to the strange quark mass as determined by the 0. We find 
My = 1.003(20) (40)GeV (second error comes from the uncertainty in scale which we take 
from the string tension) compared to the experimental = 1.019GeV. Hence our light 
quark masses span roughly the range from 0.7m s to 1.3m s . Clearly one would ideally like 
to go to much smaller quark masses. For our lightest mass we encountered one exceptional 
configuration in an ensemble of initially 200 configurations. This stopped us from attempting 
to further decrease the quark mass and reduced the number of usable configurations to 199. 
For the heavy quark mass we employed one value a s M = 4, tuned to be close to the b quark 
mass using the B s meson mass as experimental input. Fig.l shows a t Mki n for the B s meson 
extracted from B correlators with different spatial momenta. 

Mkm ~ 2 5E( P ) (10) 

with 5E(p) = E B (p) — E B (0). E B denotes the falloff energy of B correlators and differs from 
the total energy of the B meson since the NRQCD action does not include a rest mass term. 
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This distinction is irrelevant for the difference Eb{p) — Eb{0). Data from different momenta 
all give results consistent with the experimental value for Mb s and also with perturbative 
expectations. The latter is based on 

M%? = Z m M -E + E B (0). (11) 

The full horizontal line in Fig.l gives the one- loop result for a t M^ '. The two dotted lines are 
estimates of errors due to higher order corrections. In ref. [[| it was found that semileptonic 
form factors in the b quark region are not very sensitive to the heavy quark mass. Hence 
we believe errors coming from inadequate tuning of a s Mo are small compared to our other 
systematic errors. 

The "speed of light" coefficient C in the light quark action, eq.(f|), was fixed using pion 
dispersion relations. We found it sufficient to consider the lowest pion momentum at one 
value of the bare mass (a t m = 0.033). This fixed Cq = 0.94 which works well for all 5 light 
quark masses. Fig. 2 shows the speed of light C(p) defined as, 



CM = ^£>_iM (12) 

for three different light quark masses. One sees that the relativistic dispersion relation 
holds well, within errors, for momenta up to at least 1.5GeV. In our analysis of three-point 
correlators and extraction of form factors we will use continuum relativistic formulas for 
E n (p). In this work it was fairly painless to fix Co nonperturbatively. Had we used one- 
loop perturbation theory ||T5| we would not have been too far off with Cq~ 1oop = 0.91(3). 



Once two-loop results are known, we will probably be able to dispense with nonperturbative 
tunings of parameters such as Cq or the renormalized anisotropy x an d anisotropic actions 
will be just as easy to handle as isotropic ones [|18[ . 



For completeness we show chiral extrapolations for m p and Eg(0) in Fig. 3. This leads to the 
(l/a s ) p quoted above and to a B s — Bd splitting of 86(13)MeV. The experimental B s — Bd 
mass difference is 90.2(2.2)MeV. 



IV. ANALYSIS OF THREE-POINT CORRELATORS 



The first step in semileptonic B decay simulations on the lattice is to calculate the three- 
point correlator, 

Cf\p B ^t B) t) = XXe-^e'C^-^OlM^aO Vf(t,y) ^(0)|0). (13) 

x y 

$t and & B are interpolating operators used to create the pion or B meson respectively. 

is the dimensionless Euclidean space lattice heavy-light vector current. The continuum 
Minkowski space is related to V„ via 

V, = afZ Vfi £(ai) V^. (14) 
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£ (/i) is the conversion factor between Minkowski and Euclidean space 7-matrices and Zy is 
the heavy-light current matching coefficient. Perturbative estimates of Zy^ are discussed in 
an appendix. 

We work exclusively with B mesons decaying at rest (j>b = 0). t B is fixed at ts = 23 for 
time-forward runs and ts = 25 in time- reversed runs. 

We have found it convenient to rescale as follows, 

CJ?W) = m Cf\p B = 0,p*,t B ,t) e E B^ (15) 

where Eg is the ground state B meson falloff energy, obtained from fits to two-point 
correlators (in a bootstrap analysis E% is obtained separately for each bootstrap ensemble). 
Cjp is then fit to, 

£ E A jie - ESh e- E B^ I e 4 x) (*B-t) (16) 



-Ei l h 



1 

N' 

+ Y^A 2l e- E - )t e-^ E B- E B^B-t) (17) 
1 



+ 



We use constrained (Bayesian) fits of ref. to fit a single correlator to the multi-exponential 
form on the RHS of eq. flT7|) . Details are given in the appendix. Fits are carried out with £ m j„ 
fixed at tmin = 1 and for various t max < £g. The number of exponentials is increased until a 
good fit to the data is obtained. For a "good fit" we generally require x 2 /d.o.f. < 1. In a few 
cases we accept x 2 /d.o.f as high as ~ 1.2 but even in those instances the Q- values are larger 
than 0.2 . Sample fit results for (Vq) are shown in Figs. 4 and 5. For pion momenta (001), 
(Oil) and (111) and using t max = 17 ~ 21 good fits were obtained with just the first line on 
the RHS of ([17]) and > 2. The figures show results for N w = 3 and t max = 19. For the 
zero momentum case an additional excited B exponential from the second line in (|T7j ) was 
needed for good fits with t max between 17 ~ 21. Fig.4 plots the case A^- = 3 and N' w = 1. 
The fact that we can fit the data for t max only slightly below ts (the source point for the 
B meson) with just one or two Eg exponentials tells us that excited B states are highly 
suppressed in the three-point correlators. We believe this is a physical effect indicating that 
form factors for semileptonic decays of excited B mesons (proportional to Aji in ( |TTD with 
j > 1) are suppressed relative to ground state form factors. Figs. 6 and 7 show sample fits 
to matrix elements of the spatial component of the vector current, (V^). Again one can go 
out to t max ~ 21 with Nb = 1 for most momenta and with iV^ = 2 for momentum (001). 

Our goal is to extract the amplitude An, the groundstate contribution to the Cjj*\ From 
An one can then determine the B meson semileptonic decay form factors (see next sec- 
tion). One consistency check on An is to verify that its associated exponential factors, 
-e v t p -e b (t fl -t)^ i nvo ive the correct groundstate energies. For Eb this is put in by hand 



7 



through our rescaling in ( |I5|) and the ansatz of (|1T[) . Since we are always dealing with zero 
momentum B mesons for which can be determined accurately from two-point corre- 
lators, we believe this is a sensible way to proceed. For E n , for which we need results for 
various momenta, one possibility is to do simultaneous fits to two- and three-point correla- 



tors ensuring that the same set of energies appear in both correlators [19|j. We have opted 
not to force the E^p in the two correlators to be equal in this way, but to do separate fits 
and use consistency between the two independent extractions of E^\p) as a check on our 
fits, on our fit ansatz (|TTD, and on our choices for t max , Nb, etc. Fig. 8 shows groundstate 
pion energies E v (p) extracted from two-point and either (V L ) or (V fc L ) three-point correla- 
tors. One sees good agreement between the different determinations. One also notices that 
for higher momenta three-point correlators provide more accurate energies than two-point 
correlators. Of course, whether this happens or not depends on the smearings used in the 
pion interpolating operator. 

In order to include Kqcd/M corrections to the heavy-light currents we have looked at the 
matrix element of 

v ^ L = m-~^^ Q (18) 

(see ref. (2T)] for a complete list of 1/M current corrections). Figs. 9 and 10 show sample fits, 



similar to Figs. 4 and 5, for (V L ). We find that the groundstate contribution from Vq L 
is only a small fraction (1.5 -2.5%) of the leading order contribution from V L . Fig. 11 shows 
the ratio between An(Vo ) and Au(V L ). We superimpose the 0(a s /a s M) power law 
correction that must be subtracted from the (Vq ) matrix element to obtain the physical 



0(Aqcd/M) relativistic correction | 2l| . This is given by the full horizontal line, the dotted 
lines representing our estimate of uncertainties in the a s /a s M power law subtractions. We 
see that the matrix element is consistent with being 100% power law. The 1/M corrections 
to the spatial component of the vector current, (V^" ), is found to be even smaller (at the 
1% level) relative to the leading order (V fc L ). Since (V^) must be proportional to the pion 
momentum p n ^ (for pg = 0), whereas is sensitive mainly to the 6-quark momentum 
inside the initial B meson at rest, one would expect (V^' L ) to be very small. We could not 
come up with a similar plausibility argument as to why the temporal component, (Vq 1 ^' 1 ), 
should be small as well. Because the tree-level 1/M current matrix elements are so small 
and of the same order of magnitude as 0(a s /a s M) power law corrections, we have opted 
not to include them in our final analysis. Uncalculated 0(a 2 s / a s M) corrections could easily 
switch the sign of their contributions. We are dropping terms that are 1-3% of the leading 
order contributions, effects that are much smaller than the 0{a 2 s ) systematic errors we will 
be assigning to the present calculation. 

Our results for (V^' L ) disagree with those in ref. || where a much larger, 10 - 20%, contribu- 
tion from V^'' L is reported. The origin of this discrepancy is not understood at the present 
time. Nevertheless, in the next section we will see that our final results for form factors agree 
very well with those of ref. [0. The global fits used there to carry out chiral extrapolations 
must be compensating in part for the differences in the 1/M current corrections. 
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V. RESULTS FOR FORM FACTORS 



The groundstate amplitudes Au(V^) extracted from fits to three-point correlators in the 
previous section are related to the continuum matrix element of interest as 



(irM\W\B(p B = 0)) 



2 x /E n M B Z Vu 



£!p anc l are fixed from ix — ix and B — B correlators. 



Ee-^(0|$.(t,f)$t(0)|0) = E^ 



(0 



+ e 



-E«(T-t) 



(19) 



(20) 



E(o|$ B (t,f)$ko)|o) = E^V^ 



(21) 



The standard form factors f+(q 2 ) and fo(q 2 ) are defined through (q^ = p B — p%) 
(Tt(p«)\V»\B(j> B )) = U{q 2 
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M% -ml 

2 $ 



+ / (g V"* -, ,,„ 



(22) 



Following Fermilab [Q] we have found it convenient to introduce other form factors /n and 
defined as 



(vr(^)l V» \B(p B )) = y/m^yfr +plh 



with 



Pb 



M b 



(1,0) 



(23) 

(24) 
(25) 



E^ is the pion energy in the B meson rest frame and the last expressions in ( p4"D and ( |25| ) 
are similarly the four vectors and p^_ in this frame. The form factors f» and f± are useful 
since (again in the B rest frame) they are simply related to the three-point correlators 
for fi = and fi — k respectively. One has 



ft 



MV L ) 



2E^Z 



IT ^Vo 



(26) 



and 



fi 



A ^M1 J2E~ Z Vk / p w k 



(27) 



Once /|| and f± are determined, /+ and fo can then be obtained from, 
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1 



1 



(M B - E w ) fx 



(28) 



y/2MB~ 



fo 



[(M B -E„)f l{ + (Et 



Oh] 



(29) 



(M 2 - ml) 



From these formulas one sees that / + will be dominated by f±, i.e. by the matrix element 
of Vfc, and fo by f\\ or the matrix element of Vq. 

In Fig. 12 we show the form factors /+ and fo with the light quark mass fixed at 
a t m = 0.033, a value which is close to the strange quark mass. The pion momenta span 
(000), (001), (Oil), (111), (002) and (112) in units of 2n/L s a s . One sees that statistical errors 
are reasonable down to about q 2 = 16GeV 2 . More work is required if one wants to go further 
away from the zero recoil point. 

Figs. 13 and 14 show chiral extrapolations at fixed pion momentum. We have tried linear and 
constant fits to either all 5 data points or to just the last 3 points. The full and dotted lines 
in Figs. 13 and 14 give some idea of the spread in fit results. These differences are included in 
the chiral extrapolation systematic errors that we quote. With the present statistics it is not 
sensible to try more sophisticated fits. Much smaller statistical errors and data at smaller 
light quark masses are required to search for chiral log or square root type behavior. We also 
believe it is premature to try fitting to HQET and/or chiral perturbation theory inspired 
model ansaetze. In their plot of f» and f± versus am q the Fermilab collaboration, working 
at smaller quark masses than in this article, finds an upward curvature as one decreases the 
light quark mass JE]. We cannot rule out or verify such behavior with our present data. Fig. 
15 gives form factors for the physical case B — > n Iv. One sees that errors have increased 
significantly over those in Fig. 12. Furthermore we now include pion momenta only up to 
(111). Larger momenta lead to chiral extrapolation errors that are too large to make such 
data points meaningful. 

In Fig. 16 we compare our results to those by other lattice groups [fi]-||J|. One sees that agree- 
ment for f+{q 2 ) is good among all collaborations. For fo{q 2 ) we agree best with the JLQCD 
collaboration [^] and are slightly below the results of remaining groups. The Fermilab [JTJ] 
and JLQCD collaborations @ are the two other groups that simulate directly at the b quark 
mass, so it is worthwhile making further comparisons with their work. We do so for the two 
form factors (f\ + f'2) and fi used by JLQCD, which are closely related to f\\ and f±. 



Figs. 17 and 18 show comparisons between the three collaborations. The form factors are 
plotted as a function of E n , the relation between the two variables q 2 and E n being q 2 = 
Mg + m 2 — 2MsE n . The Fermilab results for (fi + f 2 ) are considerably higher than those 
from the other two collaborations. The main reason for this difference seems to come from 
the upward curvature, mentioned above, that Fermilab sees in their plots such as Figs. 13 
and 14 of form factors versus the light quark (or the pion) mass. Neither JLQCD nor the 
present work has sufficient accuracy at low enough quark masses to see this trend and more 



( /l + /2 )= * /, 




(30) 
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calculations are required to resolve this issue. One should note that soft pion theorems, 
valid in the limit m n — > and p n — ► 0, would dictate 

[A + / 2 ]|^o = ^-v / ^- ( 31 ) 

The higher Fermilab results in Fig. 17 are consistent with this relation while JLQCD's and 
our results are too low. 



VI. SUMMARY 

We have studied semileptonic B meson decays using highly improved gauge and quark 
actions on anisotropic lattices. We developed constrained fitting methods for analysing 
three-point correlators and extracting groundstate amplitudes in a controlled way. Our final 
results for form factors agree with previous lattice results. 

Our data points in Figs. 16, 17 and 18 include the main systematic errors. Allowing for 
8% discretization, 4% relativistic, 8% higher order perturbative and 2% mass tuning cor- 
rections, we estimate ~12% systematic errors from all sources other than quenching and 
chiral extrapolation. This is to be compared with the 10 — 15% chiral extrapolation errors 
already shown in Fig. 15. One realizes that accurate semileptonic form factor results will only 
be attainable if uncertainties coming from chiral extrapolations are brought under control. 
To overcome this obstacle, we have initiated a program to study heavy-light physics with 
improved staggered light quarks ||11|| . Simulations can be carried out with much smaller 



quark masses using this light quark action. The experience acquired in the present work 
and the analysis techniques that have been developed for three-point correlators will play an 
important role there. For instance, with staggered light quarks two-point and three-point 
correlators have time oscillating contributions which must be taken into account in fits. The 
only way to obtain groundstate contributions to three-point correlators will be through fit- 
ting them directly, as was done in the present article. Taking ratios of three- and two-point 
correlators will be of no use in simulations with staggered light quarks. Other theoreti- 
cal developments, such as better understanding of chiral perturbation theory for staggered 
fermions |22]| and the use of "Moving NRQCD" [|12j should further aid accurate semileptonic 



form factor determinations in the future. 
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APPENDIX A: CONSTRAINED FITTING 



In this appendix we give some details of our constrained fits to two- and three-point corre- 
lators. The general approach, within the context of lattice simulations, is described in ref. 
||. In lattice simulations one typically starts with numerical data for some correlator G(t), 
averaged over configurations, which one wants to fit to a theoretical expectation G t h{t) to 
extract energies, amplitudes or matrix elements. Examples of G t h(t) would be the RHS's of 
(|T7|), or (^T|), which we can generically write as, 

G th (t) = 1 £A n e- E » t . (Al) 

n 

Denoting the fit parameters A n and E n collectively as acj, one has 

G th (t) = G th (t J {a J }). (A2) 

Conventional fits are carried out by minimizing the x 2 , 

x\{^}) = £[G(t) - G th (t, {«,})] a'} [Gj¥) - G th (t', {a,})] (A3) 
t,t' 

with respect to the parameters {ctj}. cr^ 1 is the inverse of the correlation matrix, 



a t , t , = G(t)G(t>) - G(t) G(t>) (A4) 

Depending on the quality of the data, only a few low lying energies and amplitudes will be 
constrained by the data. If one includes too many terms in (|A1|) , the unconstrained E^s 
and A n 's for higher n can wander all over the place and start to destabilize the fits. 

"Constrained fits" were proposed in ref. || to get around this problem. One augments the 
conventional \ 2 with a term, x 2 vrior , which prevents fit parameters that are not constrained 
by the data from taking on "unreasonable" unphysical values. 



2 2 — 2 _i_ 2 ( kK\ 

X * Xaug X ' Xyriori \ 



with 



2 (q! ' — ^_/') 2 

Xprior = 5^ ' (^6) 

J a 3 

In this scheme each parameter ctj has its set of "priors", a.j and &j, and x\ U g i s designed to 
favor a.j values in the range otj ± Oy The replacement of eq . (|A5|) can be justified within the 
framework of Bayesian statistics and Bayes' theorem. It implies using a Gaussian apriori 
distribution for the parameters {<X/}. For parameters otj that are determined by the data, 
adding x 2 W ior nas minimal effect on the final fit value, as long as <5j is not made too small. In 
the present work we always set &j = ctj. One is dealing with very wide Gaussians and hence 
very unrestrictive priors. We have checked that changing a to 75% or 50% or even 25% of cej 
does not change results for data-determined fit parameters. Choices for the central values 
ctj are made based on preliminary fits and physics input about typical level splittings in the 
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system under study. Again, if the a/s are wide enough final results for data-determined 
parameters are not sensitive to precise values of the <x,'s. 

The method is best illustrated by an explicit example. Table II lists priors used in fits to pion 
two-point correlators for our next to lightest quark mass. Figs. 19 and 20 show fit results for 
the groundstate energy as a function of the number of cosh's for several pion momenta. The 
numbers below the data points show the y 2 /d.o.f. for the fits. The fit range is shown on the 
top left corner of the plots. One sees that good fits are obtained for N cos h > 4. The fancy 
stars show bootstrap fit results (bootstrap methods within constrained fits are discussed 
below). The priors E\ and A\ and E 2 were chosen from preliminary fits or by looking at 
effective mass plots. We have checked that changing them by factors of 2 moves fit results for 
groundstate energies and amplitudes by much less than their fit errors, (precise definition 
of fit errors will be given below when we discuss bootstrap methods and bootstrap errors). 
In Fig. 21 we show what happens if one changes the priors for the higher states (n > 2) from 
[E n = E 1 + (n-l)x 0.3] to [Ex + (n - 1) x 0.2], or from A n = 0.05 to A n = 0.08. 

One sees from Figs. 19 and 20 that once sufficient number of exponentials (cosh's) are in- 
cluded, fit results stabilize. We then fix N cos h and carry out bootstrap fits for our final 
analysis. For instance, for pion two-point correlators we choose N cos h = 4. 

In a bootstrap analysis involving constrained fits one first creates a certain number (we 
choose nboot=200) of bootstrap ensembles in the usual way. For each bootstrap ensemble a 
different prior value &j is picked at random for each j according to a Gaussian distribution 
about a central value c^o with width 5j. In bootstrap fits, Table II should be viewed as 
giving values for {<x/,o} rather than for {ctj} and we set Oj = ctj^. Fits are carried out for 
each of the nboot bootstrap ensembles using x 2 aug with the {ctj} for that ensemble. In order 
to get a bootstrap average and bootstrap errors one sorts the nboot fit values according 
to size and discards the top and bottom 16%. We take the average of the remaining 68% 
as our bootstrap average and one half of the difference between the largest and smallest 
values within the 68% as our bootstrap error. The fancy squares in Figs. 19 and 20 give 
bootstrap results calculated this way. One sees very good agreement between bootstrap and 
non-bootstrap single fits. For the latter, errors are calculated from the square root of the 
diagonal elements of the covariance matrix C, defined through 

Given the very different definition of errors and the fact that in the bootstrap fits very 
different priors are being used compared to in the single fits (where {ctj} = {ctj t o}), we find 
the consistency between the two types of fits very reassuring. In Figs. 22 and 23 we show fit 
results for the groundstate amplitudes, An, contributing to the (Vo) threepoint correlator 
(see eqns.(|T7|) and ( |i~9"D for definition of An). The N n = 3 fits are those that went into the 
plots of Figs. 4 and 5. One again sees good agreement between single and bootstrap fits. 
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APPENDIX B: ONE-LOOP PERTURB ATIVE MATCHING 



In this appendix we summarize the perturbative calculations necessary to match the 
NRQCD/D234 heavy-light vector current to its continuum QCD counterpart at the 0(a s ) 
and 0(a s /a s M) level. The formalism is described in detail in ref. [p0[| . We have generalized 
those calculations to include anisotropic lattices, improved glue and a more highly improved 
light quark action. The D234 one-loop self energy corrections have already been calculated 
in ref. [[TJJ for these more complicated lattices and glue actions. In the present work we do 
not include 0(a s a s ) or 0(a s Aq ^j D ) terms in the action or in the currents. Hence, in the 
notation of ref. only the Coo and Cio elements of the mixing matrix are required, in 

addition to the heavy quark self energy. The relation between current matrix elements (V^) 
in continuum QCD and the matrix elements (V^) evaluated on the lattice, is given to this 
order by, 



1 



z® 



{[l + a s p J (Vf) + (V^ L ) sub } 



(Bl) 



with 



Po„ 



B„ 



■ (C q + C Q )-C 



00,^ 



(B2) 



and 



(V^ L ) sub = (V^ L ) - a s ( 10 ,» (V, L ). 



(B3) 



The second term proportional to Cio,n * s the 0(a s /a s M) power law subtraction term plotted 
in Fig. 11. C q and Cq are the one- loop light and heavy quark wave function renormalizations, 
is the tree-level light quark wave function renormalization, 



Z (0) 



(a t m) 2 + 2(a t m)x + 1 



(B4) 



and is given by, 



B = - 

7T 



-- + ln(a s M) 



Bh = — 



7T 



-— + ln(a s M) 



(B5) 



In Table III we list one-loop results for C q , Cq, £oo,/i and Cio,^ 



Cq, Coo.aj an d Cio,At f° r a sM = 4.0 and massless light 
quarks. We work in general gauge (a g = 1 and a g = correspond to Feynman and Landau 
gauges respectively) and use gauge invariance of £io !A( and the combination [|(C , 9 +C'q)+Coo,m] 
as checks on our calculations. In Table III we present only the IR finite parts of C q , Cq and 
Coo,£t- The IR divergent pieces cancel between the lattice and continuum parts of the matching 
calculation. The Landau gauge results have smaller numerical integration errors since both 
the light quark wave function renormalization and the heavy-light vertex correction are IR 
finite in this gauge. 

For reasons described in the text, we do not include (V^) su b in our final results. The 
matching factors Zy^ of eqs.(|T9|), (|26|) and ( p7|) are then given by, 



14 



z v» = -1=^ [1 + a s p J (B6) 
V Zq 

We have used a s ~ 0.25(5) in our perturbative matching. These values are close to av(2/a s ) 
estimated on isotropic lattices with unimproved glue. Systematic errors assigned to higher 
order perturbative corrections should cover this uncertainty in a s . 
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TABLES 



TABLE I. Simulation Details. 



lattice size 12 3 x 48 

# configs 199 

P 2.4 

Landau link u n s =0.7868 u t =0.9771 

Xo 3.0 

X = a s /a t 2.71(3) 

C 0.94 

aj 1 1.20(5) GeV 

a t m 0.023 - 0.043 

P/V 0.62 - 0.76 

a s M 4.0 



TABLE II. Priors used in pion correlators for atm = 0.028 









n = 1 






mom = (000) 


0.20 


0.07 


(001) 


0.28 


0.04 


(011) 


0.34 


0.03 


(111) 


0.40 


0.02 


(002) 


0.50 


0.02 


(112) 


0.54 


0.02 


n > 1 


E 1 + (n-l)x 0.3 


0.05 
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TABLE III. One-loop perturbative coefficients for a s M = 4.0, massless light quarks and 
X = 2.71. The coefficients Cq, C g , Coo,/*) Cio,/* arm Po,/* are defined in Appendix B. a 9 is the 
gauge parameter. 





^0 






«<? = 1 


a g = 


a g = 1 


a 9 = 


Cq 

Cq 


-0.066(3) 


0.435(1) 






Coo,/* 


0.629(1) 


0.1285 


0.506(1) 


0.0070 


Cl0,/t 


-0.096 


-0.096 


0.054 


0.054 


[5(^0 + ^) 

+ Coo,/*] 


0.606(3) 


0.606(2) 


0.483(3) 


0.485(2) 


PO./i 


-0.244(3) 


-0.244(2) 


-0.334(3) 


-0.332(2) 
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FIGURES 
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FIG. 1. atMkin derived from correlators with different momenta. The full horizontal line gives 
the one-loop perturbative estimate. The two dotted horizontal lines indicate perturbative errors. 
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FIG. 2. C(p) versus the pion momentum for three light quark masses. 
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FIG. 3. Chiral extrapolation of the p mass and of Eb(0). For atm p both linear (full line) and 
quadratic (dotted line) extrapolations are shown. 
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FIG. 4. Fits to 3-point correlators. All fits have 3 exponentials coming in from the left and 1 or 
2 exponentials from the right. Both the fit and the data have been multiplied by e En t e E B (*s-t) f or 
presentation purposes. Results are shown for strange type light quarks, i.e. for the third (middle) 
light quark mass out of a total of 5. In the upper right corner we show the momentum of the pion. 
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FIG. 6. same as Fig.4 for Vfc 
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FIG. 7. same as Fig.4 for 
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FIG. 8. Consistency test for pion energies extracted from different correlators 
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FIG. 11. Ratio of groundstate contributions to (Vq 1 ^' L ) and (V^) for several pion momenta. 
The horizontal line shows the one-loop 0(a s /a s M) power law subtraction term for the (Vq 1 ^' L ) 
matrix element. 
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FIG. 12. The form factors /+ and /o for the light quark mass fixed at the strange quark mass. 
Only statistical errors are shown. 
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FIG. 13. Chiral extrapolations of the form factor /y at fixed pion momentum. Constant and 
linear fits were carried out to all 5 or to the last 3 data points. The full and dotted lines give some 
idea of spread in fit results. 
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FIG. 14. Same as Fig. 13 for the form factor 
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FIG. 15. The form factors /+ and /o after chiral extrapolation to the physical pion. Statistical 
and chiral extrapolation errors shown. 
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FIG. 16. Comparison with other groups. Statistical, chiral extrapolation and other systematic 
errors included. To avoid too much clutter we do not include errors for the Fermilab and UKQCD 
data points. They are comparable to those of other groups. 
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FIG. 17. Comparison with the Fermilab and JLQCD collaborations for the form factor (/1+/2). 
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FIG. 18. Comparison with the Fermilab and JLQCD collaborations for the form factor /2. 
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FIG. 19. Fit results c for h groundstate pion energies versus N cos h. Pion momentum is shown on 

upper right and the fit range on upper left corners. The fancy star shows bootstrap fit results. The 

numbers below data points give \ 2 /d.o.f. 
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FIG. 20. Same as Fig. 19 for pion momenta (111) and (112). 



34 



Pion fits (1-23) 



0.21 



(000) 

o Table II 

O <5E n>2 = 0.2 



0.20 — 



0.5 



w 0.4 

c0 



0.3 



Pion fits (1-23) 



(111) 

o Table II 

O <5E n>3 = 0.2 - 

□ A n>2 = 0.08 



]_g I i i i i I i i i i I i i i i I i i i l 

02468 0246 

FIG. 21. Compariso c n°s s of fits using priors of Table II with fits after chang^mg 1 the n > 2 priors 
as indicated. 
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FIG. 22. Fit results for the groundstate amplitude An from the (Vo) threepoint correlator. 
The numbers below the data points give x 2 /d.o.f . 
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FIG. 23. Same as Fig.22 for momenta (011) and (111). 
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